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("^ The key ingredients of systematic QCD analysis of MPI are discussed. 

,^ 1 Introduction 

. This contribution presents the QCD approach to studying MPI that I am involved in together 

^'^ with Boris Blok, Lonya Frankfurt and Mark Strikman. The basic new ideas of our approach 

1—^ have been announced in a short article [1]. A detailed pQCD analysis of the main contributions 

i-Q to double hard parton collisions, as well as evolution of the emerging generalized double parton 

*7^ distributions, 2GPD's, can be found in 2 . 

^ I An excess of jets + photon events in the back-to-back kinematics have signaled the presence 

^ of double parton collisions in the Tevatron experiments IHIS] . This phenomenon cannot be 

explained within a naive independent parton approximation. A model of two-proton correlations 
inside the proton that is capable of explaining the magnitude of the MPI contribution will be 
presented in [B]. 

Our approach to the problem of MPI is in certain sense opposite to that developed by 
Tevatron experiments. The CDF group in the pioneering study [3] has formulated a noble 
r~^ quest of extracting the MPI contribution without referring to either QCD theory or even to 

f^ MC event generators. Such a puristic approach has successfully performed a noble task of 

fy-s establishing the presence double hard collisions. However, it does not help much in developing 

f^ the theory of the phenomenon that would stem from the first QCD principles at the "hard" 

OA end, and exploit known phenomenology of hadron interactions at the "soft" end. Meantime, 

^"^ such a theory is necessary for quantifying two-parton correlations in the proton and for better 

^ understanding of the underlying physics of collider experiments. 

Following the pioneering work of Refs. [71 [5], a number of theoretical papers on multipar- 
^ ton interactions appeared in recent years [HI HOI IHl [HI [13] ■ They were based on the parton 

C^ model and geometrical picture of collisions in the impact parameter space. This topic is be- 

ing intensively discussed in view of the LHC program [Ml [15] ; Monte Carlo event generators 
that produce multiple parton collisions are being developed [16| , I17 | fT8 ] . In our view, however, 
important elements of QCD that are necessary for theoretical understanding of the multiple 
hard interactions issue have not been properly taken into account by above-mentioned intuitive 
approaches. 

More recently, theoretical papers exploring the nature and properties of double parton dis- 
tributions and discussing their QCD evolution have appeared [l][T9 l [20 l l21j. 

The problem of theoretical approach to MPI is, sort of, educational: both the probabilistic 
picture, the MC generator technology is based upon, and even the Feynman diagram tech- 
nique, when used in the momentum space, prove to be inadequate for careful analysis and 
understanding of the physics of multiple collisions. 
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2 Hidden reefs of MPI analysis 



A careful approach to MPI phenomena uncovers a number of unconventional features. 

Thus, in order to be able to trace the relative distance between the partons, one has to use 
the mixed longitudinal momentum-impact parameter representation which, in the momentum 
language, reduces to introduction of a mismatch between the transverse momentum of the 
parton in the amplitude and that of the same parton in the amplitude conjugated. 

Another unusual feature of the multiple collision analysis may look confusing at the first 
sight. It is the fact that, even at the tree level, the amplitude describing the double hard 
interaction contains additional integrations over longitudinal momentum components; more 
precisely — over the difference of the (large) light-cone momentum components of the two 



partons originating from the same incident hadron (see Section 2.2 1. 



2.1 Transverse structure 

Cross section of two-parton collision can be cast in the following intuitively clear form: 

daC*) = j d?pij d?p2 f{xi,pi)f{x2,P2) ■ da^P^; B = p^ - pi. 



(1) 



Here vectors p are the transverse positions of incident hadrons 1 and 2 with respect to the point 
where the two partons interact; their difference is the impact parameter of the hadron collision 
B. Local parton density / is the square of the wave function: f{xi,pi) = ip{xi, pi) ■(/'^(s^i, Pi), 
or, in the momentum representation, 



fixuPi) 



d^h^ 



■i>{xi,ki±_) 



2U 



d'k 



2/j\xi,k[^^)-e'P'^^'^-^'^^l 



ii,/,t 



(27r2) ^^^' '^' J (27r2) 
Substituting into (II]) and integrating over transverse coordinates gives 



(2) 



d^k^ 



■4'ix,k±_) 



^ ^_l,/,tr^ J,'. W / W2^Xki-kl) = 



(M- 



'ilj'{x,k'j_) X / d^pc 



'^'^^ V'(x,k^)Vnx,k^). (3) 



(27r)2 



This shows that one could have written an 
answer for the cross section in terms of mo- 
menta from the start, treating incident objects 
as plane waves with given (longitudinal and 
transverse) momentum components. 

The situation is different when two pairs of 
partons collide. Indeed, in this case the trans- 
verse coordinates of four partons are not inde- 
pendent but are related, see Fig. [T] 

Pi - P3 = P2 - P4 = S. 




Figure 1: Geometry of four-parton collision 



To see how does this condition affect the momentum picture, we introduce inclusive two- 
parton probability density distribution, and turn to the momentum representation in analogy 
with Eq. ^. 
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1=1 



12=3" fe=3 i=l '^) 

■ i^{xi,Pi;X2,P2]X3,p3; . . . Xn, Pn) i'Hxi: Pi; X2, P2; X3, ps; .. .Xn,Pn): 



Then, integrations over (pi + P2), {ps + pa) and (pi — P2) + (ps — P4) produce, respectively, 
A:i±-/c'i_L = -(fc2±-fc2_L) = A, A:3_L-fc_3_L = -(fc4±-fc4±) = A and A = -A. (5) 

The presence of the relation 

Pi- Pz- {P2 - P4) = 

leaves the transverse momcntuni parameter A 
arbitrary. 

We conclude that, in the language of mo- 
menta, a mismatch appears between transverse 
momenta of the parton in the wave function 
and the wave function conjugated. This mis- 
match is the same for all four participating par- 
tons as shown in Fig. [2j 

Thus, the new variable A is an intrinsic 
part of the two- parton correlation function in Figure 2: Shifts in parton momenta 

momentum space, that was dubbed in [1] "the generalized double parton distribution", 2GPD: 

[2]Dfixi,x2;qlqlA). (6) 

Here a, 6 mark parton species, h — the hadron, and qf, (j| (within the usual logic of parton 
distributions) stand for the corresponding "hardness scales" : the upper limits of (logarithmic) 
integrations over parton transverse momenta, kfj_ and fc|j^ . This is what concerns the transverse 
space structure. 




2.2 Longitudinal structure 

Now we shall look at longitudinal momenta 
of participating partons. 

To elucidate the problem that one encoun- 
ters here it is instructive to examine the case 
when a parton "0" from one hadron virtually 
splits into two, "1", "2", which offspring par- 
tons enter two hard interactions with partons 
"3" and "4" from the second hadron. 

This situation is shown in Fig. [Sj where 
black blobs mark two hard interactions that 




Figure 3: On-mass-shell singularity 



produce some large mass final state systems (intermediate bosons, pairs of large transverse 
momentum jets, etc.). 
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Fig. [3] is a tree amplitude. This means that knowing the momenta of incident partons 0, 
3, 4, and the 4-monienta of the produced final state systems Q13 and Q24:, one unambiguously 
determines the momenta of the virtual state partons 1 and 2. 

The problem is, at certain values of longitudinal momenta of incident partons one of the 
intermediate partons can go on mass shell. For example, the parton line "2" in Fig. [31 The 
amplitude develops a strange singularity right inside the physical region of the external momenta 
(fco,fc|,fc4 < 0, (3i3 and (5^4 positive). 

This singularity has been discussed in the literature more than once (see, e.g., p2 l I23 j ). 
However, to the best of my knowledge, its meaning and significance remained unclear before 
the explanation that we gave in [T]. 

Let us go back to the usual hard process picture, e.g., to DIS scattering. Consider perturba- 
tive splitting of an incident parton "0" into "1" and "2" , of which the former experiences hard 
scattering (gets hit by a lepton with a large momentum transfer Q^), while the latter goes into 
the final state. 



Figure 4: "On-mass-shell" partons in deep inelastic scattering 

When calculating the DIS cross section (see Fig. [4]) we put the parton "2" on mass shell 
and do not trace its fate. It may split developing a final state jet, it may propagate as a "real 
particle" at macroscopically large distances (confinement does not concern us here) and might 
eventually enter another hard interaction. This is exactly what happens with the diagram of 
Fig. [3] and where lies an explanation of the origin of that disturbing singularity. 

What happens is the following. The singularity appears at definite momenta of incident 
particles and, in particular, of partons "3" and "4" . Definite momenta mean plain waves. 
But plain waves are not localized in space-time, so that the distance between the two hard 
interactions in Fig. [3] is not known and can be in fact arbitrarily large. If the hard scattering of 
"1" and "3" would occur in the LHC tunnel, and the collision of "2" and "4" — in, say. Gran 
Sasso, then the presence of the singularity is natural; it would correspond to free propagation 
of the particle "2" between Geneva and L'Aquila. Such a scenario is possible, but this is not 
what we are looking for: we intend instead to study the situation when "3" and "4" belong to 
the same proton! 

To assure spatial localization inside one hadron, one has to construct a wave packet by smear- 
ing over the relative longitudinal momentum of the two partons. Importantly, the kinematics of 
the process determines only the sum of the light-cone momentum components {(3^ + Pi)- So one 
is allowed — and has to — introduce an integration over the difference of the two momenta, 
/33 — /34, at the amplitude level. This smearing eliminates the singularity of the diagram Fig. [3] 
the integral reduces to the residue at the pole of the propagator "2" . 

4 MPI@LHC2011 



3 Generalized double parton distributions 

We have chosen to study production of two pairs of large transverse momentum jets as an 
example of double hard interaction. The corresponding cross section is conveniently represented 
as a product of cross sections of two independent collisions normalized by the factor S that has 
dimension of area: 

da^^^ _da{xi,X2) da{x3,Xi) 1 1 _ J j^ D{xi,X2; A)D{x3,X4;-A) 
dhdh~ dfi df2 ^ S'' S~ D{xi)D{x2)D{x3)D{x4) ' *- ^ 

The quantity S is often referred to in the literature as "effective cross section" . However, a cross 
section, by definition, depends on interaction strength while S characterizes transverse area of 
two-parton correlation in a hadron and longitudinal correlation between the partons (see |6j) 

« 7SI -^ = 1^27^ = 7^ -1^27^ ^ -JT^J^- -527^- (») 



dtidt2 \Q^J R^Q^ Q6 J12Q2 ^t^^f^ J12Q2 

So, this cross section turns out to be power suppressed as compared with that of the 2 — >■ 4 jet 
production mechanism when one gets four jets out of 2-parton collision at an expense of two 
additional QCD emissions, da'^^^^K 

3.1 4- and 3-parton collisons 

The value of the correlation radius is determined by convergence properties of the A integral 
in (It]). If two partons are taken directly from the non-perturbative hadron wave function, 
2GPD=[2]£'(a;i, X2), the correlation area is of the order of the transverse size of the hadron. 

There is an additional contribution to the MPI cross section due to collision of three partons. 
In this case the numerator of (l7| has a mixed structure: 

li]D{xi,x2;A)[2]D{x3,x4;-A) + i2]D{xi,x2; A)ii]D{x3,X4;-A). 

Here [ij-D stands for 2-parton distribution involving perturbative parton splitting, as in the 
upper part of Fig [31 Being a small-distance correlation, mi I? depends on A only logarithmically 
(via parton evolution effects). As a result, the A integration in (l7| becomes broader. In spite 
of this geometrical enhancement, the 3 — >■ 4 contribution turns out to be numerically small at 
Tevatron energies {x > 10^'^) but may become significant at the LHC [6j. 

Evolution equations that incorporate into 2GPDS all-order radiative QCD effects in the 
leading coUinear approximation are described in [2]. 

3.2 A four-parton or a two-parton collision? 

There were discussions in the literature whether the process of double parton splitting shown 
in Fig. [5] should be looked upon as a 4-parton collision. 

On one hand, it looks indeed as two hard interactions of four partons. On the other hand, 
such a diagram naturally appears as a loop correction to a "normal" 2 — > 4 QCD process when 
one goes beyond the tree approximation. The question of potential double counting was raised. 
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Figure 5: Loop in a double hard process 



The process displayed in Fig. [5] is a 
"product" of two small-distance correlations, 
[^D{xi,X2)ii]D{xs,X4). Since [ijl? practically 
does not depend on A, the integral in (It]) for- 
mally diverges. 

This means that this double hard inter- 
action is not a MPI, in our interpretation 
{ [3 H] ) • It lacks a characteristic feature of MPI, 
namely a power enhancement of the differen- 
tial cross section in the back-to-back kinemat- 
ics, Qi3,(924 < Q^ (see below). 

This is a loop correction that belongs to 
2-^4 background and has to be subtracted in 
a search for MPI. 



3.3 Modeling 2GPD 

The first natural step is an approximation of independent partons, which allows one to relate 
2GPD with known objects, namely 

D{xi,X2, ql ql, A) ~ Gix,,ql;A)G{x2, ql,A). (9) 

Here G is the non-forward parton correlator (known as generalized parton distribution, GPD) 
that determines, e.g., hard vector meson production at HERA (Fig. [6]). 

The GPD, on its turn, can be modeled as 
G{xi,qlA) ~ Dixi,ql) x F2giA^), (10) 



with D — usual one-parton distribution deter- 
mining DIS structure functions and F — the 
two-gluon form factor of the hadron. 

The latter is a non-perturbative object; it 
falls fast with the "momentum transfer" A^. 




x'i,ki + A 



Figure 6: GPD in VM electroproduction 



This form factor can be parametrized differently. For example, by a dipole formula: 

i2\ 



^2,(A^) = 1 



(11) 



Here nil is an effective parameter whose value extracted from HERA data lies in the ballpark 
of TO2(a;_0.03,Q2^3GeV2)~ 1.1 GeV^. 

A simplistic approximation of independent partons does not answer the call: it fails to 
explain a factor 2 enhancement of back-to-back 4-jet production observed by Tevatron experi- 
ments [2113]. So, intra-hadron correlations between partons have to be taken into account. One 
does not know much about them a priori. However, certain information about non-perturbative 
2-parton correlations can be extracted from phenomenology of inelastic diffraction, which allows 
one to construct a viable model for the 2 GPD of a nucleon ^. 
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4 Differential distribution in back-to-back kinematics 

Four-parton interaction is a "higher twist" eventuahty. The fact that the total MPI cross section 
is power suppressed as compared with the 2 — >■ 4 cross section does not mean that 4 — >■ 4 and 
3 — )■ 4 colhsions are impossible to access at high Q^ . 

There is an essential difference between the two 4-jet production mechanisms. Namely, in 
2 — >■ 4 processes the final jets form a "hedgehog" , while double parton collisions, 4 — ?► 4, produce 
two pairs of nearly back-to-back jets. Actually, in the back-to-back kinematics the two channels 
become comparable. 





Figure 7: Enhanced kinematical configurations of four jets 

Differential distributions due to 4 — > 4 and 3 — > 4 processes exhibit double collinear en- 
hancement: they peak at small jet imbalances, Sff. ^ pfj_ ~ p|j_, Fig. [Zl 



da 



dtidt2 d'^di3d'^S2i I dt\dt2 
da I da 



da a 

oc 



2 c2 ' ^ik = PiJ- +Pk±; 



sfss: 



24 



dtidt2 d^(5i3(i2(524 / dtidt2 



gas2 



, d^ = ,5?3 ~ Sl^ » (<5i3 + <524)' 



(12a) 



(12b) 



Structure of singularities displayed in (12a) — independent enhancements in two pair imbal- 
ances — is typical for 4 — !> 4 processes. 






Figure 8: Origin of 6'^ singularity in Eq. (12b I 
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The "end-point" contributions due to 3 -^> 4 configurations with no QCD emissions between 



the parton sphtting (0 — >■ 1 + 2 in Fig.^ and the two hard vertices is enhanced as (12b). 



Singularities in Eqs. (12) get smeared by double logarithmic Sudakov form factors of the 



partons involved, depending on ratios of proper scales, see [2]- 

5 Conclusions 

QCD approach to MPI leads to the notion of generalized double parton distributions, 2GPDS. 
Higher order logarithmic QCD corrections to 2GPDS can be assembled via parton evolution 
equations derived in [U |2] in the leading coUinear approximation. Detailed formulae for total 
cross sections and differential distributions of four jet production in the back-to-back kinematics 
can be found in [2]. 

In order to reliably extract MPI contributions and get hold of parton correlations inside 
nucleon, one has to use a different strategy from that developed and promoted by Tevatron 
experiments |S1 El H] . Tevatron methods were based on measurement of angular correlations 
between jet imbalance momenta proposed in [21]. Such characteristics, however, are sensitive 
to non-perturbative physics and are strongly affected by experimental efficiencies of jet recon- 
struction. They are difficult (if at all possible) to control theoretically and should be replaced 
by studies of correlations in transverse momenta rather than angles. 
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